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Abstract. We consider the Dirichlet problem for the focusing 
NLS equation on the half-line, with given Schwartz initial data 
and boundary data q(0,t) equal to an exponentially decaying per- 
turbation u(t) of the periodic boundary data ae 2lut+le at x = 0. 
It is known from PDE theory that this is a well-posed problem 
(for fixed initial data and fixed u). On the other hand, the associ- 
ated inverse scattering transform formalism involves the Neumann 
boundary value for x = 0. Thus the implementation of this for- 
malism requires the understanding of the " Dirichlet-to-Neumann" 
map which characterises the associated Neumann boundary value. 

We consider this map in an indirect way: we postulate a certain 
Riemann-Hilbert problem and then prove that the solution of the 
initial-boundary value problem for the focusing NLS constructed 
through this Riemann-Hilbert problem satisfies all the required 
properties: the data q(x,0) are Schwartz and q(0,t) — ae 2luJt+le is 
exponentially decaying. 

More specifically, we focus on the case —3a 2 < oj < a 2 . By 
considering a large class of appropriate scattering data for the t- 
problem, we provide solutions of the above Dirichlet problem such 
that the data q x (0, t) is given by an exponentially decaying pertur- 
bation of the function 2iabe 2lut+l,i , where uj = a 2 — 2b 2 , b > 0. 

On the other hand for periodic data exactly equal to a e 2luJt+le 
at x = 0, in the case %■ < oj, the data q x (0,t) is given (exactly) 

by the different function 2abe 2liAjt+lt , where u = ^ + 2b 2 , b > 0. 
In other words, the Dirichlet to Neumann map is unstable in the 
sense that exponentially decaying perturbations of the boundary 
data q(0,t) can lead to completely different data q x (0,t). 
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1. Introduction 

We are interested in classical solutions of the following initial-boundary 
value problem 



iq t (x,t) + q xx (x,t) + 2\q(x,t)\ 2 q(x,t) = 0, x > 0, t > 0, 
(1.1) q(x, 0) = qo(x), < x < oo, 

q(Q,t) = g {t), < t < oo, 



where the function go( a; ) belongs to the Schwartz class and g (t) = 
ae 2%u}t+ie _|_ u (t^ w here a > 0, uj, e are real, u(t) decays exponentially as 
t — ► oo, and the compatibility condition qo(0) = go(0) is satisfied. We 
will assume here that —3a 2 < uj < a 2 . 

It is known [2] that there exists a unique classical solution of this 
problem (for fixed u,qo). On the other hand, the inverse scattering 
transform formalism developed in ([B], [7], pQ), in addition to qo(x) 
and g (t) also requires the function gi(t) = q x (Q, t) for < t < oo. The 
general methodology of [5] is applied to the problem (1.1) in pQ, where 
it is assumed that the unknown function is the sum of 2iabe 2luJt+le 
(where uo = a 2 — 2b 2 , b > 0) and a Schwartz function. (In fact, [1] 
consider only the case u = 0, but their results go through even if tt is 
exponentially decaying, or, say, Schwartz.) 

It is known [8] that this assumption is not always true. Here is a 
counterexample for a 2 / 2 < u> < a 2 , when u = 0: consider the exact 
one breather solution 



;i.2) q(x,t) = 2ne l 



Airi 2 t 



cosh2n(x — xq) 



Clearly q(x, 0) is Schwartz and g(0, t) = ae 2tut+le where a = cos ^ VXo ^ 
and uj = 2n 2 . So uj > a 2 /2 but q x {0,t) = 2abe 2iujt+ie where b = 
n tanh{2nxo) and b = uj/2 — a 2 /4. 

The aim of this paper is to prove that this assumption is correct for at 
least some exponentially decaying u. Since the above counterexample 
shows that it is not true for all such u we deduce that the Dirichlet to 
Neumann map for the above initial boundary value problem (the map 
that takes q(0,i) to q x (0,t) ) is highly unstable. 
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2. A RlEMANN-HlLBERT PROBLEM 



The focusing NLS equation admits the Lax pair 



(2.1a) 
(2.1b) 



H x + ik[a 3 , fi] = Q(x,t)n, 
fi t + 2ik 2 [a 3 , jj] = Q(x, t, k)/j,, 



where a 3 
(2.2) 

Q(x,t) -- 



diag(l, 



1), 





q(x,t) 



q(x,t) 




Q(x, £, k) = 2kQ - iQ x a 3 + i\q\ 2 a 3 . 



A novel method for analysing initial boundary value problems for in- 
tegrable nonlinear PDEs was introduced in [5]. This method, which 
is based on the simultaneous spectral analysis of both the x-problem 
and the t-problem in the Lax pair, was rigorously implemented to the 
NLS on the half-line with Schwartz initial and boundary conditions in 
[7J. In the problem (1.1) the initial data are of Schwartz class, thus the 
scattering and inverse scattering of the x-problem is classical and goes 
back to the original investigations of Gelfand, Levitan and Marchenko 
(see [JJ). On the other hand, the boundary values at x = are per- 
turbations of finite-zone functions, thus the spectral analysis of the 
t-problem involves aspects of the finite-zone theory. In this paper we 
will consider the simplest possible case of zero-zone data. 

The zero- zone solution of NLS, namely q(x, t) = q p (x, t) = ae 2lbx + 2lujt + 
gives rise to the Dirichlet data ae 2lu)t+le and also yields q x (0,t) = 
2iabe 2iuJt+it . 

Now, let b be defined by u = a 2 — 2b 2 , b > 0. We will assume here 
that a 2 — u > and b 2 < 2a 2 . Let Q(k) be the function defined as 

(2.3) Q(k) = 2{k - b)X(k), X{k) = ^{k + bf + a 2 . 

Following [Ij we consider the two-sheeted Riemann surface X defined 
by the function Q(k). Our Riemann-Hilbert problem will be defined on 
X. We also consider the oriented contour E defined by Imfl(k) = 0, 
see Figure 1. (This is Figure 9 of [IJ with some contours reoriented.) 

One easily sees that the curve S consists of two copies of the real 
line and an analytic arc r U T connecting the two branch points E = 
—b + ia,E = b — ia and the two infinities ooi and oo 2 (on the two sheets 



of X). 



4 



S. KAMVISSIS AND A.S.FOKAS 




Figure 1. The two-sheeted Riemann surface X. 

£ defines a partition of the sphere X into D\, D 2 , D 3 , D 4 , where 
D x = {Imk > 0,ImQ(k) > 0}, 
D 2 = {Imk > 0, ImQ(k) < 0}, 
D 3 = {Imk < 0, ImQ(k) > 0}, 
D 4 = {Imk < 0, ImQ(k) < 0}. 



Next, define the following matrices 




H(t, k) = exp(iuxT3t)E(k)exp(—iuHT3t), 
k) = H(t, k)exp{i[uj - Q(k)]a 3 t). 



Let the functions a(k) and b{k) be the (classical) scattering data for 
the function qo{x) defined in [7j. All we need to know here is that a(k) 
is smooth for k real and can be analytically extended in the upper half- 
plane, with a(k) = 1 + 0(l/k) as k —>■ oo. Similarly, b(k) is a Schwartz 
function for k real which can be extended to the upper half-plane such 
that b(k) = 0(1 /k) as k —>■ oo. Furthermore, \a 2 \ + \b 2 \ = 1 for k real 
and a can have at most a finite number of simple zeros in the complex 
plane, say ki, k 2 , , k n , with Im(kj) > 0, j = 1, ...,n. 

Let the functions A, B be functions satisfying the following condi- 
tions: 

(i) The functions A(k), B(k) are analytic in D\\JD 3) bounded in D\U 
D 3 and satisfy the following asymptotics A(k) = 1 + 0(1/ k), B(k) = 
0(l/k) as k — > oo. 
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(ii) b(k)A(k) - a{k)B{k) = in D x . This is the so-called global 
relation. 

(hi) A{k)A{k) + B{k)B{k) = 1, A(k) ^ 0, fceE. 
We will now dehne a Riemann-Hilbert problem in X, with jump data 
given in terms of a, b, A, B, following [lj. 
We dehne the matrices 

a(k) b(k)\ 
-b(k) a{k) J ' 



(2.6) 



s(k) 



(2.7) S(k) = 

and G(k) = s~ l (k)S(k). Let 



A(k)_ B(k)\ 
-B(k) A(k)J 



p{k) = 
r(k) 



Gn(kY 
b(k) 



a(k) ' 
c(k) = p{k) — r{k). 

Consider now the following Riemann-Hilbert problem with the jump 
contour S: 

M_(x, t, k) = M+(x, t, k)J(x, t,k), fceS, 
lim k ^ OQl M{x,t, k) = I, 



[2.8) 



where 
(2.9) 

J(x, t, k) 
J(x, t, k) 



1 

r ^ e 2i(kx+(n(k)-uj)t) 
1 

J(x, t, k) 



_^(J^ e -2i(kx+(n{k)-u))t) 

l + \r(k))\ 2 

_pffc\ e -2i{kx+(n(k)-u>)t) 

l + \p(k))\ 2 

1 

c + Q^ e 2i{kx+(U(k)-uj)t) y 

1 

_ c -^ e 2i(kx+(n(k)-Lu)t) ^ 

I _c+n c \ e -2i(kx+{n(k)-u)t) 
1 



1 





nA e -2i(kx+(p.(k)-u>)t) 

1 



k e R upper , 

k e R lower , 
k e r 12 , 
k e r 21 , 
k e f 12 , 
k e f 21 . 



Here c + and c_ are boundary values of the function c which is analytic 
in D 2 . 
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res k = kj [M(x,t, k)]i = im)e 2i{ - kx+ ^-^[M{x, t, kj)} 2 , kj E D u 

res k=Zj [M{x,t,k)]i = im 2 e 2 < kx+ ^-^\M{x,t, Zj )} 2 , z 3 E D 2 , 

res k= - Zj [M(x,t,k)] 2 = -imJe-^+W^-^tMfo t, %)]i, z 3 E D 3 , 

res k=k .[M(x,t,k)} 2 = -im)e 2 < lx+ ^-^[M{x, t, kj)] 1} k A E D,, 



Theorem 2.1. The above Riemann-Hilbert problem admits a unique 
solution. 

The theorem follows immediately from the so-called vanishing lemma 
extended to the surface X [9] by employing the symmetries of the jump 
J. Although the vanishing lemma applies to holomorphic Riemann- 
Hilbert problems, the above meromorphic Riemann-Hilbert problem 
can be easily transformed to a holomorphic Riemann-Hilbert problem 
as in [3] by adding small loops around the poles and changing variables 
inside the loops (see also [6], [7]). 

3. Asymptotic analysis of the Riemann-Hilbert problem 

The analysis in section 3.3 of [T] shows that the Riemann-Hilbert 
problem above gives rise to a solution of the focusing NLS in the first 
quadrant. Furthermore the initial data q(x, 0) are equal to q because 
of the definition of a,b. What is not a priori clear is that q(Q, t) = 
g (t)+u(t) andg x (0,t) = 2iabe 2iuJt+ie +v(t), where u, v are exponentially 
decaying at infinity. 

This is the main result of this paper. 

Theorem 3.1. Define q(x,t) = 2ilim k _> OQ1 kMi 2 (x, t, k) where M 12 is 
the (12) entry of the solution of the above Riemann-Hilbert problem. 
Then q(x,t) solves the focusing NLS equation in the first quadrant, 



withq(x,0) = q Q (x), q(0,t) = g (t) + u(t), q x (0,t) = 2iabe 2iujt+ie +v(t) 



where u(t),v(t) are exponentially decaying at infinity. 

PROOF: Follows from the asymptotic analysis of the Riemann-Hilbert 
problem (for data a,b,A,B), as t — ► oo. From section 3.3 of [1] we 
have that the Riemann-Hilbert problem above reduces to the following 
Riemann-Hilbert problem when x = 0: 



where 



(2.11) 



m) 



(ib( k j)-^r( k j)) X > m ) = -res k=Zj c(k), 



m) 
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(3.1) 
where 



M®(t, k) = M { *\t, k)J {t \t, k), fceE, 
lirrik^MVfak) =1, 

I B(k) c -2i(n~u))t s 



(3-2) J^\t,k)=[ m : m ^ )t A{k ) , l.JbGE, 

\ A(k)A(k) J 

where the superscript + denotes the limit from the +side of the contour 
and the superscript — denotes the limit from the —side of the contour. 

The following asymptotic analysis will show that as t — > oo, we 
recover the pure zero-zone solution. 

Theorem 3.2. Up to an exponentially small error, the Riemann-Hilbert 
problem for is asymptotically (as t — > oo) equivalent to the trivial 
Riemann-Hilbert problem which has no jump. By this we mean that 
lim k ^ oc {k n 'Mi2 ) is exponentially small for n = 1,2,.... 

Proof. Note the factorization of J w on R upper U f : 

J®(t, k) = J up J l °, 
l/D_ B(k)A(k)D_e- 2i ^-^ 
(3.3) ~V D ~ 

D+ 
A(k)B(k)(D + )- 1 e 2 ^ n -^ t 1/D, 

and where D solves the scalar problem 

D + = D_A(k)A(k), k e R upper U f 

and satisfies limk^ ool D{k) = 1. This factorization follows from the 
identity A(k)A(k) + B(k)B{k) = 1 for k e E. 

Similarly, note the factorization of on R lower U T: 

J {t \t,k) = G lo G up , 



where J up 
and J l ° 



(3.4) 



1 

where G ° = \ B(k) r2 m- U )t i 

A(k) e 

1 g(fc) p -2i(Q-a;)r 
£<«P _ [ - 1 - A(fc) c 

' 1 



For the asymptotic analysis we must deform our Riemann-Hilbert prob- 
lem in small lenses with boundaries consisting of the different compo- 
nents of R U T and slight deformations of these components. 
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For example we consider the oriented contours C 1,up and C 1,l ° from 
ooi to 002 on the upper sheet of the Riemann surface slightly deforming 
the real line, with C 1,up lying in Di and C l,l ° lying in Z) 4 , and denote 
the corresponding lenses D l up and Di i D in a way that dDi up = C 1,up U 
M. upper and dD 2}Up = C 2 ' up U M<w er . We construct similar lenses around 

r r ]$} ower 

We define O as follows: 

0(t, k) = M®(t, k)J up , k e Du , 

(3 ' 5) 0(t,k) = M^(t,k)(J lo )-\ keD liUp . 

Similarly for the other lenses. Note that O is piecewise analytic off E 
only if A, B are analytic in the appropriate lenses. This is not assumed 
to be generality true, but A, B can always be approximated by analytic 
functions in a way that the overall error due to the substitution of A, B 
by their analytic approximations is exponentially small as t — > 00 (see 

We now observe that the off-diagonal entries of the jump matrix 
for O are uniformly exponentially small. On the other hand, the 
diagonal entries are uniformly bounded. So, according to standard 
asymptotic analysis of Riemann- Hilbert factorization problems [I], it 
follows that, up to an exponentially small error, O is given by the 
solution of a problem with diagonal jump, which in turn reduces to the 
scalar problem for D. The off-diagonal entries of thus have to be 
exponentially small in t, to all orders in k. □ 

The limiting Riemann-Hilbert problem is trivial and corresponds to 
the purely zero-zone solution of NLS. Using the formulae q(0,t) = 
Um^JpikMn^t)] and fe (0,t) = Um k ^ oo [4k 2 M 12 (0,t)+2iq(0,t)kM 22 (0,t)} 
we see that u(t),v(t) are actually exponentially small. Theorem 3.1 is 
thus proved. In fact, using similar formulae for J^-f(t), i = 1,2,3, .... 
in terms of M'*' it is possible to show that u, v are Schwartz functions. 
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